Electromagnetic wave scattering by Schwarzschild black holes 
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We analyze the scattering of a planar monochromatic electromagnetic wave incident upon a 
Schwarzschild black hole. We obtain accurate numerical results from the partial wave method 
for the electromagnetic scattering cross section, and show that they are in excellent agreement with 
analytical approximations. The scattering of electromagnetic waves is compared with the scattering 
of scalar, spinor and gravitational waves. We present a unified picture of the scattering of all massless 
fields for the first time. 
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Black holes are thought to be efficient catalysts for the 
liberation of rest-mass energy. As such, black holes are 
implicated in the most energetic phenomena in the known 
universe (e.g. gamma ray bursts). On the other hand, 
after a turbulent youth, many black holes settle into a 
quiescent old age. Some estimates suggest there may be 
up to a billion quiescent stellar-mass black holes within 
our galaxy [T]. Their existence may be inferred from, 
for example, the transient lensing of background sources; 
a handful of events have so far been observed [5]. A 
possibility for future consideration is that quiescent black 
holes may be indirectly identified from the 'fingerprint' 
they leave on radiation that impinges upon them. 

Over the last four decades, some clues about the prop- 
erties of any such 'fingerprint' have been uncovered. For 
example, a time-dependent perturbation incident upon 
a black hole will excite characteristic damped ringing in 
response. The frequencies and decay rates of the ringing 
are linked to the well-studied quasinormal mode spec- 
trum [3]. Black holes illuminated by long-lasting planar 
radiation will create interference patterns, and rotating 
black holes will create distinctive polarization patterns 
[4]. Both effects depend strongly on the ratio of horizon 
size to wavelength. Hence, it is conceivable that future 
gravitational-wave detectors may be able to identify the 
fingerprint from rapid and distinctive variations across a 
narrow frequency band. Nevertheless, inferring the pres- 
ence of quiescent black holes from such clues must remain 
a challenge for future decades. 

Scattering by black holes is of foundational inter- 
est in both black hole physics [Sj and scattering the- 
ory [S]. Many authors have studied the simplest time- 
independent scenario, in which a black hole is subject to 
a long-lasting, monochromatic beam of radiation. Here, 
the key dimensionless quantity is the ratio r/j/A where 
is the horizon size of the black hole, and A is the wave- 



length of the incident wave. The interference pattern 
depends also on the spin s of the perturbing field, with 
s = 0, 1/2, 1 and 2 corresponding to scalar, neutrino, 
electromagnetic and gravitational fields, respectively. 

To the best of our knowledge, the first paper outlin- 
ing a calculation of wave scattering cross section in the 
spacetime of a black hole was published by Matzner [7J in 
the late sixties. Since then, planar wave scattering from 
black holes has received much attention, especially in 
Schwarzschild and Kerr spacetimes (see Refs. [SI [HI IS] for 
comprehensive accounts on the subject). Let us briefly 
review a sample of the literature for the simplest case, 
the Schwarzschild black hole, for which the scattering of 
monochromatic fields of all spins (s = 0, 1/2, 1 and 2) has 
been studied through the years. The case of scalar waves 
(s = 0) was extensively studied by Sanchez [THIIIT], both 
analytically and numerically, and an accurate numerical 
study was later performed by Andersson [12. . Fermion 
(s = 1/2) scattering by a Schwarzschild black hole was 
the subject of a recent study [1] , in which the authors also 
elucidated the effect of non-zero field mass. The case of 
electromagnetic waves (s = 1) was studied analytically 
by Mashoon [T3| and Fabbri [iT , and some results were 
obtained in the low- and high-frequency limits. Gravita- 
tional waves (s = 2) were the first to be studied in black 
hole scattering [IS], and are the subject of old [5J[T^ and 
new [HI [H] works. 

In this letter we present the first detailed numerical in- 
vestigation of the scattering of an electromagnetic plane 
wave by a Schwarzschild black hole. This work fills a 
gap in the literature, and complements recent numerical 
studies of the scalar [T^ , fcrmionic ^ , and gravitational 
[T5] cases. We take this opportunity to present a unified 
picture of all four fields, for the first time. 

We use natural units with c = G = \ and the metric 
signature (H ). 
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The line element of Schwarzschild spacetime can be 
written as 



for X :$> Vs {r ^ Ts). Here h\^\x) denote the spherical 
Bessel functions of the third kind |23], and and 
T^jP are the reflexion and transmission coefficients, re- 



ds ^ f{r)dt - [/(r)] dr -r [dd + sin ed(j) ), (1) spectively, which satisfy 
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where /(r) = 1 — 2M/r, with M being the black hole 
mass. The Schwarzschild solution describes static and 
chargeless black holes, with event horizon at = 2M . 

The Lagrangian density of the electromagnetic field in 
the modified Feynman gauge is 120] 



with g = det(5^^), G = V^'A^ + K^A^ and = 
(0, df /dr, 0, 0) . The equations of motion are found to be 



(2) 



The two physical polarizations in Schwarzschild space- 
time can be written as 
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(3) 



(4) 



with w > 0, and I ^ 1. (See, e. g., Ref. [20] for a dis- 
cussion on the possible solutions of Eq. In Eqs. ([s]) 

and Yim and F^f™ (a = 6*, 0) are the scalar and vector 
spherical harmonics ^21j , respectively, and ip^^ (r) satisfy 
the following equation 

{u:^ - V{r)) ^ti (r) + (^f^^ti (r) ) - 0, (5) 
with A = /, //, where the effective potential is V{r) = 

mi+i)/r% 

To evaluate the solutions of the Eq. ([s]) in the asymp- 
totical limits we use the Wheeler coordinate, defined as 
X = r + rsln {r/r^ — 1), and rewrite Eq. ([s]) as 



{u'-V)vti{^) + ^Hii^))^0- 



(6) 



For the computation of the scattering cross section we 
need only to consider modes incoming from the past null 
infinity J^^ . For these modes, the asymptotic solutions 
of Eq. ^ are 122] 



for 



ftii^) 
— oo (r — !■ Ts) and 



A A rjiX —iuJX 
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(7) 



A' 



(8) 



normalization constant which is not important for the 
scattering properties. 

The phase shifts are related to the reflexion coefficient 

by 



(9) 



For Schwarzschild black holes, the phase shifts of the 
two different physical polarizations are the same, i. e. 

The differential electromagnetic scattering cross sec- 
tion is 1241 



da 

dn 



1 

4^ 



^ 2/+1 



rp/(cos( 



de 



Pi\cos9) 



(10) 



where P™(cos6') are the associated Legendre functions. 
Note that Eq. ( |lO| takes into account the contributions 
from the two physical polarizations, and it is valid for 
both linearly and circularly polarized waves. The polar- 
ization properties of the initial wave remain unchanged 
in the scattering by non-rotating black holes [13]. 

For small angles, this scattering cross section is the 
same for the massless scalar and electromagnetic fields, 
and it is given by [13l [14] 



da 

dn 



16M2 



(11) 



In fact, the same behavior for small angles is obtained 
for massless fermionic and gravitational fields scattered 
in Schwarzschild spacetime [H [5j [17] . 

The glory approximation for scattering of electromag- 
netic waves by a Schwarzschild black hole can be deter- 
mined using the strong field approximation for the de- 
fiection angle (which was first obtained by Darwin [25] ) 
together with the general glory formula [26], namely 



da 

dn 



2Trujbi 



[J2s{ujbg sin( 



(12) 



where b is the impact parameter of the incident particle, 
Ji{x) are the Bessel functions of first kind, bg is the im- 
pact parameter for which the scattering angle is tt, and s 
is the particle spin. For the electromagnetic field {s — 1), 
the glory scattering cross section is given by 



1 da 

Jpdn 



30.75Mw[J2(5.36Ma;sin6l)]^ (13) 



where the coefficients (30.75 and 5.36) were obtained by 
solving the geodesic equation numerically |27| . 
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Schwarzschild scattering cross section for the electromagnetic field at Mm - 1 .0 



Schwarzschild scattering cross sections at A/o) = 4.0 
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Schwarzschild scattering cross section for the electromagnetic field at Mm - 2.0 
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Schwarzschild scattering cross section for the electromagnetic field at Mui - 3.0 
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FIG. 1: Electromagnetic scattering cross section of 
Schwarzschild black holes for different choices of Mlj. We 
compare our numerical results (solid lines) with the glory ap- 
proximation (dashed lines) given by Eq. (13 1, obtaining ex- 
cellent agreement for 6* « vr. 
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FIG. 2: Scattering cross section of Schwarzschild black holes 
for massless scalar (s — 0), electromagnetic (s = 1), gravi- 
tational (s = 2) and massless fermionic fields (s = 1/2) at 
Muj — 4.0. Note the log scale on the vertical axis of the lower 
plot. We see that, as all other non-zero spin fields, the elec- 
tromagnetic wave has a vanishing scattering cross section in 
the backward direction. 



In order to evaluate the electromagnetic scattering 
cross section numerically, we first solve Eq. ^ and match 
the solution with Eqs. (|8| and The scattering cross 
section for arbitrary frequencies and angles is obtained 
through Eq. (10). The numerical method used here is 

. We have 



analogous to the one described in Ref. 
employed an iterative method similar to that used in 
Refs. p!51 \T9\ to improve the numerical convergence of 
the partial wave series. 

In Fig.[T]we show the differential electromagnetic scat- 
tering cross section of Schwarzschild black holes com- 
puted numerically for different values of the incident wave 
frequency {Mw = 1,2,3,4). We also show the results 



for the glory scattering [given by Eq. ( 13 )] in each case. 
Our numerical results are in excellent agreement with the 
glory approximation for 9 n. 

The zero in the backward direction (Fig. [T]) is a conse- 
quence of the parallel-transport of the polarization vector 
along a geodesic. Consider an incoming geodesic ray in 
the z-direction which orbits the hole once, to return in the 
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opposite direction (9 = tt). Assume, without loss of gen- 
erality, it has an electric-field vector in the x-direction. 
If the ray orbits in the x-z plane then the vector will be 
reversed, whereas if the ray orbits in the y-z plane the 
vector remains unchanged. Hence, by integrating over 
the circular degeneracy (all orbital planes), there is per- 
fect cancellation. Similar arguments hold for other spins 

Eiiaol. 
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In Fig. [2]we plot the differential scattering cross section 
of Schwarzschild black holes for massless scalar (s = 0), 
massless spinor (s ~ 1/2), electromagnetic (s = 1) and 
gravitational (s = 2) waves. As expected, in the back- 
ward direction all non-zero spin massless fields have van- 
ishing cross section, whereas the zero-spin (scalar) mass- 
less field has a glory maximum at 9 = tt. We see that 
scalar (s = 0) and electromagnetic (s = 1) scatter- 
ing cross sections are very similar in the angular region 
A5° < 9 < 160°. All integer spin fields (s = 0, 1,2) be- 
have similarly for 45° < 6* < 120°. Bosonic (s = 0, 1, 2) 
and fermionic (s = 1/2) scattering cross sections oscillate 
in antiphase throughout almost all the angular range of 
Fig. |2] (except near 9 ^ tt). 

The regular oscillations in the cross sections of Fig. [2] 
can be understood semi-classically. They arise from the 
interference of rays passing in opposite senses around the 
hole. There is a 'path difference' between the rays passing 
through angles 9 and 277 — 9. A maximum (minimum) oc- 
curs when the path difference is an integer (half-integer) 
multiple of the wavelength A. Hence the angular width 
of the oscillations is inversely proportional to Muj. 

In summary, we have studied the scattering of 
a monochromatic planar electromagnetic wave by a 
Schwarzschild black hole. We have applied the partial 
wave method to obtain the differential scattering cross 
section numerically, for different values of the frequency 
of the incident plane wave and for different values of the 
scattering angle. We have presented graphs with accu- 
rate numerical results for massless fields of all spins, that 
is, scalar (s = 0), fermionic (s = 1/2), electromagnetic 
(s — 1) and gravitational (s = 2) fields. All non-zero spin 
massless fields have a vanishing scattering cross section 
in the backward direction {9 — tt), whereas the scatter- 
ing cross section of the massless scalar field has a local 
maximum. 
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